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ABSTRACT 

A coalescence model using the observed properties of pre-stellar condensations (PSCs) 
shows how an initially steep IMF that might be characteristic of primordial cloud frag- 
mentation can change into a Salpeter IMF or shallower IMF in a cluster of normal 
density after one dynamical time, even if the PSCs are collapsing on their own dy- 
namical time. The model suggests that top-heavy IMFs in some starburst clusters 
originate with PSC coalescence. 
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1 INTRODUCTION 

A recent study of observations of the stellar initial mass 
function (IMF) suggest there are systematic variations 
where the IMF gets flatter, or more top-heavy, in denser re- 
gions (Elmegreen 2004). This paper proposes that the high 
mass part of the IMF varies with density as a result of the 
coalescence of dense pre-stellar condensations (PSCs), such 
as those observed by Motte, Andre, & Neri (1998), Testi & 
Sargent (1998), Onishi et al. (2002) and Nomura & Miller 
(2004). These objects have densities in the range from 10 5 
cm" 3 to 10 7 cm -3 , and masses from 0.1 Mq to a few Mq, 
giving them sizes of ~ 10 4 AU. The largest PSCs may be 
self-gravitating (Johnstone et al. 2000, 2001; Motte et al. 
2001; Umemoto et al. 2002), as may those with stars (Fachi- 
hara et al. 2002). The size and spatial proximity of PSCs in 
a cloud core suggest their coalescence might be important 
for the highest masses (Elmegreen & Shadmehri 2003) . This 
process is modelled to illustrate its possible effect on the 
IMF. 



2 THE MODEL 

Consider the evolution of a distribution n(r, M, t) of PSCs 
with masses M and radii Rp at positions r inside a cloud 
having a mass M c and radius _R C . The cloud density profile 
is 



Pc(r) 



PcO 



1 + (r/Rco) 



(1) 



where the core radius is R c o and the central density is given 
by the cloud mass: 
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The PSCs are taken to have density profiles (Whitworth 
Ward-Thompson 2001) 

Ppo 



Pp (r P ) = 



(3) 



[1 + (r P /R P o) 

where central density p p o is assumed to be constant and the 
radius R p varies with condensation mass M p and position r 
in the cloud. The positional variation assumes the density at 
the edge of the condensation equals the ambient cloud den- 
sity, p p (ii p ) = p(r), so PSCs of the same mass are smaller 
near the cloud core. Thus the density contrast between the 
condensation center and the edge is 

C(r) = (1 + [r/R c0 ] 2 ) 



• P P o/PcO, 

and the PSC radius can be written as 

R P (r,M) = a(r)R pQ 

where 

R P o = (M/27rp p0 ) 1/3 ^arctan [a(r)] 
and 
a(r) = 



a(r) 



1 + a(r) : 



'1/3 



(4) 



(5) 



(6) 



(7) 



R p (r, M) is assumed to be the PSC radius at the beginning 
of its collapse. 

The average radius of the PSC used in the collision rate 
should in general be smaller than R p because of PSC col- 
lapse. Each PSC is assumed to collapse at a rate propor- 
tional to its central dynamical rate, (Gppo) 1 ^ 2 , so each PSC 

radius drops effectively to zero as e~"( Gpp0 ) * for param- 
eter v ~ 1. For a distribution of formation times r up to 
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the current time t, these assumptions give the average PSC 
radius for use in the collision cross section: 



< R p (r,M,t) >= 

J* R p (r, M)n(r, M, r) exp(-i/ v /Gp p() [t - r])dr 
J*n(r,M,r)dr 



(8) 



Now that the PSC radius is defined, the PSC collision 
rate and the evolution of the IMF can be calculated. Two 
components are assumed for the IMF at intermediate and 
high mass. One is intrinsic to the formation of PSCs and 
independent of coalescence. This could be a power law if 
the PSC formation process is scale-free, or it could be log- 
normal if there is a characteristic mass and some type of ran- 
dom multiplicative fragmentation or build-up. To illustrate 
how apparently power law IMFs can result from IMFs with 
a characteristic mass, the model considers a log-normal in- 
trinsic spectrum of PSC mass, F(log M)dlog M = f(M)dM 
'(to use the notation of Scalo 1986). This is taken to be 
normalized so that J A " max f(M)dM = 1 for minimum and 
maximum stellar masses M m in and M max . Power law models 
were also made and the results were qualitatively the same 
as those reported here so they are not discussed further. 

The time-evolution of the number, n, of stars between 
masses M and M + dM at time t and distance r from the 
center of the cloud is 



dn{r,M,t) _ ep c {r){Gp c \r}) 



1/2 



dt 



+0.5?? 



-f(M) 



f 

J M 



<M> 

n(r,Z,t)n(r,M -t,t)<r(£,M 



-r/n(r, M,t) 



/■Af max 

/ n 

J Mmin 



(r,Z,t)<T(Z,M)vd£, (9) 



where < M > is the average mass in the IMF, e and r\ are 
efficiencies, and a and v are the coalescence cross section 
and relative velocity. 

The first term is the spontaneous star formation rate. It 
is written to give a local rate ep (Gp c [r]) 1 '' 2 in units of mass 
per volume per time, where e is the fraction of the gas mass 
converted into star mass in each unit of dynamical time, 
(Gp c y 1/2 (see Elmegreen 2002). Typically e ~ 0.2-0.5 from 
observations of cluster formation efficiencies and timescales 
(Lada & Lada 2003; Elmegreen 2000). When divided by the 
average star mass in the IMF, < M >= J Mf(M)dM, 
it gives the number of stars formed per unit volume and 
time. The function f(M)dM is the probability that a star 
has the mass M between M and M + dM. Multiplica- 
tion of the count rate by f(M) gives the number of stars 
forming per unit volume and time in a dM mass interval 
around M. Note that the star formation mass per volume 
and time is J Mn(r, M,t)dM , which is ep (Gp c ) L ^ 2 again, 
since / Mf(M)dM =< M > cancels the < M > in the 
denominator of equation (9). 

In the composite IMF model, f(M) should have a steep 
slope at intermediate to high mass so the overall IMF has 
a steep slope (|T| > 1.7) in low density environments, as 
is apparently observed (Elmegreen 2004). The Miller-Scalo 
log-normal function is used here because it has this property. 

t In this paper, log refers to base 10 



Scalo (1986) cautioned that the Miller-Scalo function is not 
a good representation of the field-star IMF. He replaced the 
log-normal with a slightly steeper function at intermediate 
mass and a slightly shallower function at high mass. Here, 
the Miller-Scalo log-normal is used only as a starting point 
for the calculation, calling it a spontaneous IMF, and then 
it is evolved with coalescence to get a final IMF. For this 
log-normal, which is defined for equal intervals of log(M), 
the probability function f(M), which is defined for equal 
intervals of M, is given by: 



f(M) 



foe 



-B(log[M/M ]) 2 



M 



(10) 



where /o is the normalization factor making J f(M)dM = 1. 

The second term in equation (9) represents the collision 
rate to form a PSC of mass M per unit intervals of time, 
volume, and mass. The cross section a includes gravitational 
focusing, 



a{M x ,M y ) = n(R px + R p 



1 + 



2G (M x + M y 



2v 2 (R px + R p 



(11) 



The collision velocity v for PSCs is not well observed. It 
may be the same for all objects and comparable to the virial 
velocity in the whole cloud, or may be a function of distance 
from the center, t) oc r°, in analogy with some turbulence 
models. Here, v is assumed to equal a representative three- 
dimensional virial velocity, 



(3GM C /5R C ) 



1/2 



(12) 



The coefficient in front of this expression to correct for 
the non-uniform mass distribution in the cloud is relatively 
unimportant; most of its effect can be absorbed in the fac- 
tor 77. The third term in equation (9) is the rate of loss of 
condensations of mass M resulting from collisions with all 
other masses (see Field & Saslaw 1965). 

Equation (9) is solved as an initial value problem by 
assuming n(r, M, 0) = at t = 0. The number density of 
PSCs between masses M and M + dM starts out like f (M) 
for all r, but once a few PSCs form, collisions convert low 
and intermediate mass stars into high mass stars, changing 
the IMF. Because most of the cloud mass is at low density 
where collisions are relatively unimportant, the IMF is not 
affected much overall, but at small radii it can change a lot, 
becoming flatter at high mass. 

Figure 1 shows several solutions to equation (9). 
The plotted quantity is an integral over radius in the 
cloud, 2.3M J fl n(r, M, t)4nr 2 dr, where the multiplication 
by Mj log(e) = 2.3M converts the function n to log inter- 
vals for plotting. The cloud parameters are: M c — 10 4 Mq, 
Rco = 0.1 pc, and R c — 5 pc. The PSC central density is 
p p o = 10 s p cm" 3 for mean molecular weight p = 4 x 10~ 24 
gm. The collapse rate parameter for PSCs is v = 1. The 
efficiencies in the evolution equation are: e = 0.5 and 77 = 0, 
0.45 and 0.95. Larger 77 increases the effect of collisions and 
flattens the IMF; 77 = shows the spontaneous mode. The 
intrinsic IMF has parameters as in the other examples above: 
M min = M = 0.1 M Q and B = 1.08. The red, green and 
blue lines represent IMFs integrated out to 2, 5, and 10 times 
the cloud core radius, R c o', they illustrate how the cloud core 
IMF is shallower than the others. The total integration time 
is one dynamical time in the cloud core, td = (Gp c o) _1/ ^ 2 . 
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Figure 1. IMF model with collisions between pre-stellar conden- 
sations following the evolution equation (14) up to one dynam- 
ical time in the cloud core. Curves are labelled by the value of 
r\ (dashed curves have no coalescence: r\ = 0). Color indicates 
radius over which the IMF is integrated: red, green and blue in- 
tegrate out to 2, 5, and 10 times the cloud core radius. The IMF 
is shallower closer to the center of the cloud because coalescence 
is more important there. Fiducial slopes are shown by straight 
dashed lines with the T values indicated. The Salpeter function 
has r = —1.35. The IMF at intermediate-to-high mass is steep in 
the absence of collisions (T ~ —1.7), as might apply in low den- 
sity star-forming regions in the field, and it becomes shallower 
over time, especially in the cloud core, as coalescence makes more 
massive stars. 



Evidently coalescence in dense clusters can turn a steep 
intrinsic IMF into one that resembles a Salpeter function, 
with a slope of ~ —1.35 or shallower. The same spontaneous 
IMF without collisions would have a slope of 1.7 or steeper, 
depending on the mass range. 



do not necessarily imply a single scale free star formation 
mechanism (Elmegreen 2004). Our work in this paper con- 
firms this expectation by presenting a simple model for colli- 
sions between PSC in clusters and its possible effect on IMF. 
Of course, we did not intend to fit the IMF with a particu- 
lar coalescence model, but showing that the high mass part 
of IMF in a composite model as suggested by Elmegreen 
(2004), may change due to interactions between PCS. 
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3 SUMMARY 

We considered a cluster with nonuniform density profile and 
known mass and radius, in which PSCs can move approx- 
imately with virial velocity. Considering collisions between 
PSCs and their gravitational collapse, a differential equa- 
tion is presented and solved numerically which gives number 
density distribution of PSCs inside the cloud at each radius, 
mass interval and time. In order to control the collision and 
the collapse rates in a phenomenological way, some input 
parameters are introduced. We explored a wide range of in- 
put parameters and the results were qualitatively the same 
as those reported here. In some starburst regions a high 
fraction of high mass stars compared to low mass stars has 
been reported which means a high mass bias (e.g., Sternberg 
1998; Smith & Gallagher 2001). Our results suggest this shift 
in high mass part of IMFs in these starburst clusters may 
originates from PSCs coalescence. 

Elmegreen (2004) discussed how various physical pro- 
cesses operating in at least three mass regimes may lead 
to three components IMF. This point is demonstrated by 
reviewing the observations of IMF variations and using two 
simple models showing the observed power law distributions 
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